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Abstract We present WIYN SparsePak observations of the diffuse ionized gas (DIG) halo 
of NGC 89 1 . Preliminary results of an analysis of the halo velocity field reveal 
a clear gradient of the azimuthal velocity with z which agrees with results for 
the neutral gas. The magnitude of the gradient has been determined, using two 
independent methods, to be approximately 15 km s -1 kpc -1 . 

1. Observations and Data Reduction 

Data were obtained during the nights of 10-12 December 2004 with the 
SparsePak IFU (see Bershady et al. 2004; 2005). We used the 860 1/mm grating 
at order 2, which provides a spectral resolution A/AA ~ 4,900, and the 316 
1/mm grating at order 8, which provides a spectral resolution ~ 10,000. In 
the latter mode, the wavelength coverage includes the Ha, [Nil] AA6548,6583 
and [S II] AA67 16,6731 emission lines. The higher spectral resolution data are 
presented and discussed in this paper. The data were reduced in IRAF. Figure 
1 shows the SparsePak pointings on an Ha image of NGC 891. 

2. Azimuthal Velocities from Envelope Tracing 

To analyze the kinematic structure of the DIG halo of NGC 891, position- 
velocity (PV) diagrams have been constructed from the spectra of all 82 fibers 
in pointing 'H' (cf. Fig. 1). To enhance signal-to-noise, several fibers were 
averaged together for each major axis distance R at the highest z. Figure 1 
indicates which fibers were included in each z -range. Because a bright sky 
line interferes with several of the Ha profiles, and to increase signal-to-noise, 
PV diagrams were constructed of the sum of the [Nil] A6583 + [S II] A6716 
emission lines. Contour plots of the PV diagrams are included in Figure 3. 

Rotation curves were derived for each PV diagram using the envelope trac- 
ing method (e.g., Sofue & Rubin 2001). Because NGC 891 is nearly edge-on 
(i > 88°; S waters 1994), the edge of each line profile furthest from the sys- 
temic velocity (the "envelope") corresponds to gas at the line of nodes; this 
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Figure 1 Ha image of NGC 
891 from Rand et al. (1990), 
with the positions of the 
SparsePak fibers during the 
higher spectral resolution ob- 
servations (H) overlaid as 
black circles. Lower spec- 
tral resolution pointings (LI, 
L2, and L3), to be presented 
in a future paper, are plotted 
with gray circles. Fibers on 
the major axis are plotted in 
white for clarity. Ranges of z 
used to construct PV diagrams 
(see text) are labeled zl, z2, 
and z3. The rotation center is 
marked with a white star. The 
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velocity is thus approximately the azimuthal velocity at that R. The azimuthal 
velocity is found by (Sofue & Rubin 2001): 

v az = (v cnv - v sys )/ sin(i) - \Jal iSt + <T 2 gas , (1) 

where -u cnv is the velocity at the location of the envelope, ai nst is the velocity 
resolution of the instrument, and <T ga s is the velocity dispersion of the gas. 
Azimuthal velocity curves were constructed from the PV diagrams, and are 
shown in Figure 2. Extinction in the midplane prevented reliable velocities 
from being found for z ps 0". 



Figure 2 Azimuthal veloc- 
ity curves, derived from the 
PV diagrams shown in Figure 
3 using the envelope tracing 
method. Velocities are shown 
for 25" < 2 < 45" (zl; filled 
circles), 45" < z < 65" (z2; 
open squares), and z > 65" 
(z3; crosses), and are relative 
to v sys = 528 km s _1 . 

It is clear from Figure 2 that the derived azimuthal velocities decrease with 
height. A close inspection of the PV diagrams shows that this is not an effect of 
signal-to-noise variations with height. An average azimuthal velocity gradient 
of 0.80 km s^ 1 arcsec -1 (17 km s _1 kpc" 1 ) was calculated. Similar results 
were obtained by following this procedure for the [N II] A6583 emission line 




-160 -140 -120 -100 
R (arcsec) 



DIG Halo Kinematics in NGC 891 



3 



alone (14 km s kpc ). These values are in agreement with that measured 
by Fraternali et al. (2005) for the neutral component, 15 km s _1 kpc -1 . 

3. Azimuthal Velocities from PV Diagram Modeling 

The envelope tracing method is sensitive to the changing signal-to-noise ra- 
tio of the data with z. To account for this effect and the radial gas distribution, 
we have generated galaxy models. The Ha image presented in Figure 1 was 
used to obtain estimates of the radial density profile at each of the heights con- 
sidered in the model. Because the distribution of the DIG is not axisymmetric, 
the profiles were modified by hand to better match the shape of the data PV 
diagrams. The amplitude of the radial density profile was chosen such that the 
signal-to-noise in the model approximately matched that in the data. 

We then created model galaxies with these derived radial profiles using a 
version of the GIPSY task GALMOD, modified to allow for a vertical gradient 
in azimuthal velocity which begins at a height zq above the midplane: 

dv 

v(R, z) = v(R, z < z Q ) - — [\z\ - zq\ , (2) 

where dv/dz has units of [km s _1 arcsec -1 ], and for the models considered 
here, the major axis rotation curve is flat [v (R, z < zq) = 200 km s" 1 ]. The 
inclusion of the parameter zq was motivated by the results of Fraternali et al. 
(2005), who find that the gradient in the neutral component of NGC 891 starts 
at approximately z = 1.3 kpc (though that result may be a consequence of 
beam smearing). That value of zq was used for the models shown here, but 
we note that the appropriate value cannot be determined for our data (we lack 
data below z m 1.2 kpc). Artificial SparsePak observations were made of the 
models to create PV diagrams, and compared to the data (Figure 3). 

The model that appears to best match the data has dv/dz = 15 km s _1 
kpc -1 . A statistical analysis of difference images suggests that the best match 
occurs at dv/dz = 12 — 14 km s _1 kpc -1 , but because the model cannot 
perfectly reproduce the shape of the PV diagrams, comparing by eye is the 
most reliable way to determine the gradient. To ensure that our results are not 
influenced by an incorrect specification of the radial density profile, we also 
consider models with a flat radial density profile, and find that the same gradi- 
ent yields the best match. We conclude that the vertical gradient in azimuthal 
velocity for the DIG halo of NGC 891 is approximately 15 km s" 1 kpc -1 . 

4. Future Work 

To attempt to understand the origin of the azimuthal velocity gradient de- 
rived in this paper, the kinematic data will be compared with the results of an 
entirely ballistic model of disk-halo flow (cf. Collins et al. 2002). Further 
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Figure 3. Comparison between data (white) and model (black) PV diagrams for different 
values of dv/dz, at the indicated heights. Contour levels for data and models are 10, 15, 20, 25, 
30, 35, and 40a for 25" < z < 45" (zl); 5, 8, 11, 14, 17, and 20a for 45" < z < 65" (z2); 3, 
4.5, 6, 7.5, 9, 10.5, and 12a for z > 65" (z3). The systemic velocity is v ByB = 528 km s _1 . 



information regarding the kinematic structure will be obtained once the lower 
spectral resolution pointings are reduced and analyzed. Results from these 
efforts will be presented in a forthcoming paper. 
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